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a • 

We study the Langevin dynamics of the standard random heteropolymer model by mapping the 
problem to a supersymmetric field theory using the Martin-Siggia-Rose formalism. The resulting 
model is solved non-perturbatively employing a Gaussian variational approach. In constructing the 
■ solution, we assume that the chain is very long and impose the translational invariance which is 

expected to be present in the bulk of the globule by averaging over the center the of mass coordinate. 
In this way we derive equations of motion for the correlation and response functions C(t,t') and 
' R(t,t'). The order parameters are extracted from the asymptotic behavior of these functions. We 

find a dynamical phase diagram with frozen (glassy) and melted (ergodic) phases. In the glassy 
phase the system fails to reach equilibrium and exhibits aging of the type found in p-spin glasses. 
Within the approximations used in this study, the random heteropolymer model can be mapped to 
the problem of a manifold in a random potential with power law correlations. 
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I. INTRODUCTION 



Here, we analyze the kind of dynamical trapping that can occur in the standard model of random heteropolymer 



Disordered systems can be extremely hard to solve, as the example of spin glasses shows It took enormous 
effort to understand the physics of infinite-dimensional spin glasses, while the fate of finite-dimensional spin glasses 
is still debated. Certainly, the complexity of the spin glass energy landscape is the major obstacle one has to deal 
with, and there are other systems sharing this feature: standard examples are proteins H and manifolds in random 
\ potentials ||. 

The study of simplified random heteropolymer models may provide a useful first step toward understanding the 
physics of proteins. Here, a central question is whether the trapping of the protein in a valley of the rough energy 
landscape can hinder, or perhaps even prevent, folding into its native state. Something related to this scenario has 
I > actually been observed in some real proteins: the protein can be heated and then, upon re-cooling, misfold and never 
be able to find its native state 

g; » 

So far, in addition to numerical simulations, two analytic approaches have been used to solve such models: equilib- 
rium analysis employing the replica technique (see, e.g., refs. |]7|-|l3fl) and dynamical studies using Langevin dynamics 

g: 00. 

In the equilibrium approach, one studies the properties of Gibbs equilibrium. Even the simplest kind of random 
heteropolymer model can only be approached analytically in approximated ways. In ref. jl3| the model was analyzed 
with replica variational approximations, which predict ergodicity breaking at low temperature, giving one-step replica 
symmetry breaking (1RSB) for d > 2 and continuous replica symmetry breaking for d < 2. Thus, in 3D, below 
the freezing temperature, the ergodic components lie far apart from each other, and the same interstate overlap qo. 
Given the intrinsic one dimensional nature of the polymer, it has been argued that RSB could be an artifact of the 
variational approximation. [pgj We believe however, that even in this eventuality, implying that a single native state 
dominate the thermodynamics, RSB in the variational approximation is a signal of a complex energy landscape, which 
can lead to slow dynamics, with off-equilibrium behavior on long time scales. 

Models exhibiting 1RSB (such as the simple random heteropolymer mentioned above, the p-spin glass, or a manifold 
in a random potential) have been found to have different dynamic and static phase diagrams, with a dynamical energy 
density higher then the one found at equilibrium. This raises the intriguing possibility that, for a suitable range of 
temperatures and times, a heteropolymer might find itself dynamically trapped in a local state (as in the scenario 



1 



described above), while the equilibrium statistical mechanics might give no clue that this was happening. Such 
trapping would thus be an intrinsically non-equilibrium effect, and a dynamical theory is required to describe it. 

For the models with this feature, the solutions exhibit a breakdown of time-translation invariance (the correlation 
functions depend on the time since the system was quenched into the glassy state) and a breakdown of the fluctuation- 
dissipation relation (which is a fundamental characteristic of Gibbs equilibrium). Together, these properties of the 
dynamical glassy phase go under the name "aging" , and it is one of our goals here to examine the possibility of aging 
in heteropolymers. 

In this paper we consider the simple random heteropolymer model with Langevin dynamics (as in ref. fl6| , pT|| ) . 
The equations of motion are constructed in such a way that the Gibbs distribution is the stationary solution of the 
dynamics. This type of dynamical approach was used successfully in spin glass models. 

To derive closed equations of motion for correlation and response functions we resort to a Gaussian variational ansatz 
similar to the one used at equilibrium. The same approach has been used to study the problem of a manifold in a 



random potential, for both statics [p3|]24[ | and dynamics 25 2q|. In related dynamical work on random heteropolymer 



model |L6| and |21|] , the slightly different approach of Mode Coupling Theory is used. Our approach gives results fully 
coherent with the ones obtained there, although the detailed form of the phase diagram differs, due to the different 
nature of the approximation. 

The analysis of the variational equations indicates that, as expected from static treatments, the random heteropoly- 
mer model exhibits spontaneous breaking of ergodicity in a glassy phase. All these states are equally distant from 
each other; they have same interstate overlap (naturally, the self overlap is different). We also discuss the nature of 
the transition from the frozen (glassy) to the melted (ergodic) phase. Furthermore, we find that, within the Gaussian 
variational approximation that we employ, the random heteropolymer model can be mapped onto the the problem of 
a manifold in a random potential with power law correlations. 

The paper is organized as follows. Section II briefly describes the Langevin model. In section III a mapping to a 
supersymmetric (SUSY) field theory is made. The resulting action can be simplified by assuming a very long chain. 
This is discussed in section IV. Dynamical equations in SUSY notation, given in section VII, are obtained via the 
variational ansatz discussed in section V and VI. Also, in section VII, the connection of the random heteropolymer 
model to the problem of a manifold in a random potential will be commented upon. After disentangling the SUSY 
notation, one obtains dynamical equations for correlation and response functions (section VIII) . An analytical ansatz 
for solving these equations is introduced in section IX, and the solution is obtained in section X. Section XI discusses 
the ergodic phase, while in section XII the spin glass phase is analyzed. Technicalities needed to construct the full 
phase diagram are given in section XIII. 



II. THE MODEL 

The model is defined as follows. The Langevin dynamics is assumed to be governed by the Hamiltonian iT[ar], 

dx{s,t)/dt = -dH[x]/dx(s,t) + r)(s,t), (1) 

where x(s, t) is the position of chain bead s at time t. Beads are numbered continuously from s — to s — N. rj(s, t) 
is Gaussian noise: 

( V (s,t)r ) (s',t')} T =26(s-s')8(t-t')T (2) 

due to contact with a heat bath at temperature T. The Hamiltonian H[x] — Hq[x] + H ran d[x] contains a deterministic 
part Hq [x] and a random part H ranc i [x] . The Ho [x] is defined as 

rp r-N 

H Q [x] = - ds[(dx(s, t)/dsf + nx(s, tf] (3) 
2 Jo 

and describes the elastic properties of the chain and a confinement potential which fixes the density of the protein. 
The random part H ran< i describes heterogeneity of the interactions between the beads, 

1 f N 

H rand [x] = - / dsds'B sy V(x(s, t) - x(s', t)). (4) 



2 



B s s i is quenched Gaussian noise with variance B 2 : 



(Bl s ,) B =B 2 1 s>s>. (5) 
V(Ax) is a short-range potential, and for simplicity we take it to have a Gaussian form, as in ref. fujfl : 

V(Ax) = (-^) d / 2 e -( A -) 2 /2-. (6 ) 
2tt<7 

d is the dimensionality of the system, and a parameterizes the range of the potential. Large (small) er results in a 
long- (short-) range potential. In particular, for a — ► 0, V(Ax) — > <5(Ax), and we recover the potential used in [ pl| . 
Here and in the following Ax denotes the bead-to-bead distance: Ax = x(s, t) — x(s' , t) for a pair of beads s, s'. 

This model admits a stationary solution characterized by a Gibbs distribution. The equilibrium partition function 
for this solution is given by 

Dxe ~h Jo dsl(dx( S )/dsf+»x( S ) 2 ]-£ J" dsds'B BiB ,V(x(s)-x(s')) ,^ 

T in Eq.(|) ensures that the chain constraint and quadratic confinement are temperature- independent. (That is, the 
elasticity is purely entropic in origin.) The same convention was used in ref. |l6fl . This differs slightly from the work in 
ref pi] , where the elastic term had a factor /3 in front of it. Our choice ensures that for high temperatures the random 
heteropolymer behaves as a Gaussian random coil. Also, for very low temperatures, the random part of interaction 
with (3 in front dominates (/3 — > oo; the elastic and confinement terms become negligible). Thus, in principle, for 
(3 = oo, Z in Eq. (^) is dominated by minima of J dsds' B SiS 'V(x(s) — x(s')). Furthermore, in this limit there is 
nothing that would control the spatial spread of those minima, and (x 2 (s,t)}T diverges for very low temperatures. 
(This only happens when \i is held fixed. If it it adjusted appropriately, one can keep (x 2 (s,t))T fixed instead. In this 
paper, however, we will be concerned with finite-T phase transitions, not the low-T limit, so we will work with fixed 



III. MAPPING TO THE FIELD THEORY 

Using the standard Martin-Siggia-Rose formalism |27|, the dynamical average of any observable can be calculated 

as 



(0[x,x]) T = J DxDxDD£D£0{x,x)e- s[xS; t3, (8) 
with the following dynamical action: 

S[x,x,£,£\= J dtds -Tx{s,t) 2 +x{s,t) (j^x(*,t) + |^y ) 

- f dttfM^t) + jdtdsds^t)^^^^) (9) 

x, £, £ are auxiliary fields which appear in the formalism. To compactify the notation we introduce the superfield $: 

$(s, t u 0i, Ox) = x(s, h) + C(s, tx)Bi + Sifa, h) + 0x9x2(8, h), (10) 

where 6 and 9 are Grassmann (anti-commuting) variables. For X, X' S {9,9,9' ,9'}, {X, X'} = and J dXX = 1, 
the rest of the integrals being zero. In the following, for practical reasons, the more compact notation $(s, 1) = 
$(s, ti, Ox, Ox) will be used. Also, the integral symbol J" d9xd9xdtx will be denoted by J" dl. 

In supersymmetric (SUSY) notation Eqs. (||) and (^|) translate into ([ll|) and (|l2|): 

(0[$]) T = y* J D$0[$]e- s[ * 1 , (11) 
5[$] =S Q [$} + S rand [<Z>], (12) 
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where 



S [$] = 1/2 J dsdlds'd2<f>(s,l)K^' $(s'2), 

Srand^] = 1/2 / dldsds' B S>S ,V , 1) - $(V, 1)), 



and 



= (Jia^lfi , = T [ M - (9/9 S ) 2 ] - D[ 



(2) 



n( 2 )_ 97 1 I 

1 901001 ' ""-BOidt! dti 



20i- 



<9 2 d 



(13) 
(14) 

(15) 
(16) 



As noticed by De Dominicis Eg ] the expression in Eq.([ll|) is already normalized, so the average over the quenched 
random interactions B SiS i can be done directly on (|ll|): 



-(S [*]+Si [#]) 



((A[$]) t )b = J D<bA[>S>]e 
where exp(— Si[<&]) = (exp(— S ra nd[$]))B- The average over B s ^ s > can be done easily, leading to 

r i 2 

5i[*] = -B 2 /4 / dsds' /" dlV($(«, 1) - $(«', 1)) . 



(17) 



(18) 



The dynamical action S" = So + Si closely resembles the effective Hamiltonian obtained in the static replica approach 
of ref. fj]|§. (This rather general similarity between replica and SUSY treatments has been discussed in ref. p9L) 
Instead of summation over replica indices in we have J dl . Our expressions are not identical to those in , 
since we use a quadratic well potential instead of two- and three-body interaction terms to confine the polymer. Also, 
we use a Gaussian V(Ax) instead of S(Ax). 



IV. LONG CHAIN APPROXIMATION 

The Si part of the action can be further simplified. It can be rewritten in the form 

B 2 



Si 



with the notation 



= / dld2dxdyA^ 2 } (x,y)A[*l(x,y) 



where A^ and A^ are given by 

4]a / d8V(9>(8,l)-x)V($(8,2)-y) 



It is useful to transform exp(— Si) as 
r B 2 



exp 



-A^ * A® 



DQ 1 DQ 2 exp 



A { $(x,y)- 


—Si) as 






'B 2 r 


= exp 


16 1. 



(19) 
(20) 

(21) 
(22) 



{A (v) + A W) * (A^ + A&) - (A™ - A&) * (A^ - A&) 



B 



2 . 



(Qi *Qi + Q 2 * Qa) + Qi * {A (v) + A&) + iQ 2 * {A {v) - A^) 



(23) 
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Then, the dynamical generating functional F defined by 

e- F = [ D<S>e- s W+ M , (24) 



with J * $ = J dsdlJ(s, l)Q(s, 1), can be written as 

e~ F = J DQ 1 DQ 2 e~^ t > Ql * Ql+Q2 * Q ^ J D$e L l Ql ' Q "*\ (25) 



with L given by 

D 2 



L = —[Q 1 * (A^ + A&) + iQ 2 * (A^ - A®)] - S [$] + J * (26) 

So far everything was exact. A< v > and A<*> are both of order TV and for very long chains one can calculate integrals 
over Qi and Q 2 in (25) using a saddle point approximation. The saddle point equations read 

Q>*- = l(Am+AW) L , (27) 
QTi v - = \{A^-A^) Vl (28) 

where Q^i denotes the average with L taking Qi, Q 2 — > Q^' ,Q S 2 P ' ■ This leads to self-consistent equations for Q^' p ' 
and Q S 2 V ' ■ 

Thus, Eq.(p5[) can be approximated as 

e~ F « J ^$e- s '[*]+ J **, (29) 

with S'[9] = S([*] + 5 [*] and 

S[[$] = ?- [(A^) s , * (AW) S , - A< v > * (A«) s . - AW * (A< y )) 5 '] • (30) 

(A^^)s' and (A^^)s' have to be calculated self consistently with S': 

In the limit N — > 00 Eq.s( p9| , p0| , ^l| ) provide an exact description of random heteropolymer dynamics. 

V. VARIATIONAL ANSATZ 

To solve the model we proceed by using a variational ansatz, assuming that fields $ are approximately described 
by a Gaussian action 

S var = ^J dldsd2ds'$(s, l)G(s, 1; a', 2)~ 1 $(s', 2). (32) 

This approach has been widely used in statics. Here we apply it to a dynamic calculation. The goal is to calculate F 
given by Eq. (p9|) . Since the variational parameter G(s, 1 ; s' , 2) is the only quantity we are interested in, there is no 
need to keep the source J. It is convenient to write Eq. (pffi), with J = 0, formally as 

e- F = (e^ s '- s -\ ar e- F -, (33) 
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where 



J D<5>e~ s - , (.) var = e F - J D^(.)e- S - . (34) 

In usual statics, for problems without disorder, the variational approach is related to a maximum principle. The 
equivalent of Eq. |I3) leads to the inequality 



> e -«*'-s— )>.~ e -*,«.. (35) 



In the present dynamical problem, as well as in the static problem with replicas, unfortunately such a maximum 
principle is not known, and the variational free-energy cannot be claimed to be an upper bound to the true one. 
Despite that, the variational approach has been argued to give exact results in some limiting cases |23]j24|l , giving a 
justification for its use in general. 

The dynamical variational free-energy F^ yn = ((<S" — S var )) var + F var is given by 

Fdyn — Fdyn + F^J n + F^J n , (36) 

with 



F 



(i) 



d 



dyn 2 

Ffi = ~Tr\nG (38) 



J dsdlds 1 d2K{( G\i (37) 



F^ n = -^(A^) var *{A^) var . (39) 

Note that in calculating F^ n , the average ( )s> in Eq.([50|) is performed over the trial distribution (and therefore 
denoted ( ) var )- 

VI. CALCULATING (A (v) ) * (A i6) ) 

The term (A^) var * (A^) var is the only non-trivial term in F& yn . Before calculating it we simplify each of the 
factors in product further 

A^\x,y) « A^f\u) = 1/v [ d d RA^ 5 \u,R), (40) 



where v is the volume of the protein and a new coordinate system has been introduced: 

R=(x + y)/2,u=(x-y)/2. (41) 

R is the center of mass and u a relative distance coordinate. Thus, translational invariance is introduced by hand via 
Eq. ([4b"|) . This approximation is not necessary; the model could be solved without it. However, as we shall see later on, 
this approximation leads to dynamical equations which are identical to those for the random manifold model studied 
in |H|g. 

Changing integration variables from d d xd d y to d d Rd d u (the Jacobian is 2 d ) gives 

dld2d d xd d y{A ( $(x,y)) var {A ( i t l(x,y)) var « ^ J d\d2d d u J d d R 1 'A<$ (u, R') J d d R"A ( ^ 2 {u,R") (42) 

The integrals over R' and R" can be easily performed and one gets 
2 d 



v 



dl d2 d d u ds da' d d a d d (3 V(a) V((3) (6[2u-a + /3- $(s, 1) + 2)]) var (S [2u - $(.s, 1) + 2)]) var , (43) 



G 



which can be further written as 
2 d 



dld2d d udsds' d d ad d f3V(a)V(/3) . 
v J J {2Tr) a {2n) a 



d d p d d q 



J{p+q)2u e ip(l3- a ) e -^ P 2 B s 12e -^q 2 B 3 1 



where averages over S var have been evaluated as 

/ *p(*(«,i)-*(a,2))\ _ p -y 2 B{ 2 

\p Jvav — c , 

with Bf 2 is given by 

B',2 = ([*(«> 1) - *(«> 2)] 2 >„«r = G(s, 1; s, 1) + G(s, 2; a, 2) - 2G(s, 1; s, 2) 
Integrating Eq.(E3) first over u, and then over q and p finally gives 

-(2Tr)- d / 2 ( dld2dsds'd d ad d f3V{a)V{(3){B s 12 +B( 2 )- d l 2 c— ( " ~ 



exp - 



2(B( 2 + B[ 2 ) 



(44) 



(45) 



(46) 



(47) 



Eq.(ft7|) holds for any V(Ax). However, technically, it is of little use unless the integrals over a and j3 can be performed 
explicitly. For the Gaussian form for V(Ax) (see Eq.|J) it is possible to perform the integrals, and one gets 



-(4n)- d/2 I dldZdsds' \{B S 12 + b( 2 )/2 + a 



-d/2 



and finally 



with 



.(3) 



/ dld2dsds'V 



(B{ 2 + B{ 2 )/2 



V(z) 



B 1 



(z + a) 



-d/2 



B = 



B 2 N 

77 



(47T) 



-d/2 



Eq.s (H|3,||||) and © fully determine F dyn . 



(48) 



(49) 



(50) 



VII. EQUATIONS OF MOTION IN SUSY NOTATION 



Given the Fd yn , one can derive the equations of motion from the stationarity condition 



8G : 



(51) 



The most complicated term is g J ss , F^J n . From ([49]), it is 
d 



2N 



dMAdudvV + fl£ 4 )/2] 6 S 



(S13S23 + S14S24 - Si 3 S 2 4 - ^14<523)- 



Due to translational invariance in s, B\ 2 is independent of s. After dropping the index s Eq.(p2|) simplifies to 



5 -AH 

c' r dim ~ a U ss' 



SGl s 2 dyr ' 



■ [ d3V'(B 13 )-V'(B 12 ) 



(52) 



(53) 



The variations of F^P and FY' n are trivial. Using Eq. (|5l|) and (|36|) leads to 



,(2) 
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Kg -(G%)- 1 + 25 a 

which can be written as 



: J d3V(B 13 )-V{B 12 ) 



0, (54) 



K{G\i = 8 l2 8 ss , + 2 J d3V'(B 13 )(G s £ - G{(). (55) 
Due to translational invariance in the variable s it is useful to define following Fourier transforms 

Then Eq.([55"|) translates into 

T{n + fc 2 ) - D { P] G\ 2 = 5,2 + 2 J dZV'{B l3 ){G k 32 - G k 12 ). (57) 

Eq. (^) is identical to the one obtained in ref. [^6| for a Z?-dimensional manifold ip(uS) {to £ R D , <p G R d ) in a random 
potential V(ip(u>),u), where the correlations of the potential are described by 

(V((p, u>)V((f/, (J)) = -d5 D {u - oj')V [(ip - (p')/d] . (58) 

These equations of motion were derived using the Gaussian Variational Approximation (GVA), which is exact for the 
random manifold problem in d = oo. We expect the same behavior for the random heteropolymer. However, in this 
study we work at finite d, so the equations of motion are approximate. 

There have been several studies of random manifolds where V describes power law correlations as in Eq. (|5S|), 
employing static [p3p^j3l[|, and dynamical p0|,|32|,|33[ approaches: 



V(z) = (z + a) 1 ^/2(l- 7 ). (59) 

By comparing Eq.s (|50| ) and (|5^) one notices that V(z) is identical to V(z) (up to a proportionality factor B 2 ) if one 
identifies 7 = 1 + d/2. Accordingly, we conclude that, within the Gaussian variational approximation used in this 
study, random heteropolymer dynamics is identical to the dynamics of the manifold in a random potential with power 
law correlations. (We can not say anything rigorous outside the framework of the Gaussian variational approximation 
scheme, of course.) 

Furthermore, correlations of the random manifold potential are classified as short range for 7 > 2/(2 — D), and 
long range for 7 < 2/(2 — D) |2^j2^]. This classification of random manifolds helps to classify random heteropolymer 
model in the same way. Using 7=1 + d/2, the random heteropolymer has D = 1, and short range correlations 
for d > 2 and long range correlations for d < 2. (Again, this all makes sense only within the Gaussian variational 
approximation) . 



VIII. DISENTANGLING SUSY 

G\ 2 encodes 16 correlation functions, out of which only two, correlation and response function, are independent 
and nonzero: 

({x(s,ti)x(s',t2))) = C(s,t i;s ',t 2 ) = I ^-e^-^C^hM) (60) 

{{x(s, h)x(s', f 2 )» = R(s, t V ,s', t 2 ) = J ^L e ik ( s -^R k (h,t 2 ). (61) 

Also, by adding an external field term to the original Hamiltonian H[x] — > H[x] + J dsdtx(s,t)h(s,t) one gets 
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((x(Mi)i(s',^))) = 



(62) 



£/i(s',f 2 ) 

i.e. R(s, t\\ s', t 2 ) describes the response to an infinitesimal field applied at time t 2 and bead s' . Thus, G\ 2 reduces to 



G\ 2 — Ck{ti,t 2 ) + (6>i — 9 2 ) [0\R k {t 2 ,tx) — 6 2 R k (tx,t 2 )] , 

and, accordingly, with Gff = C(s, ti; s, tx), G 22 — C(s,t 2 ; s,t 2 ), and Eqs. ( |60l , |6l"l ) , one gets 

B 12 = B(t 1 ,t 2 ) - 2(0i - 62) [dxr(t 2 ,tx) - 6 2 r(t u t 2 )] , 



(63) 



(64) 



with 



and 



f dk 

B(t u t 2 ) = / ^[C k {t 1 ,t 1 ) + C k {t 2 ,t2)-2C k {t u t 2 )] 
r(h,t 2 ) = J -^R k (ti,t 2 ), 



(65) 



(66) 



After disentangling the equations of motion in SUSY notation (see Eq. 57) by using (63-|66|) gives 



[T(p + k 2 ) + d/dt]C k (t, f) = 2TR k (t', t) + 2 dsV [B{t, a)] R k {t' , 



+4 / dsV" [B(t, s)} r(t, s) [C k (t, t') - C k (s, t')} 



[TO* + k 2 ) + d/dt]R k (t, f) = 5(t - f ) + 4 / dsV" [B{t, a)\ r(t, s) [R k (t, t') - R k {s, t')] . 

Jo 



(67) 
(68) 



The equations of motion (67) and (|68|) are almost identical to the ones found in ref. [|25| (here D = 1, while in |2[ 
D = 0). 



IX. ANSATZ FOR C K (T,T') AND R K (T,T) 



These equations of motion are coupled integro-differential equations which in principle can be solved; the initial 
conditions are given by C k (Q,0) and we use Ito's convention R(t + e, t) — » 1 as e — > 0. It is well known that 



asymptotic solutions of such ec 
equations analytically. pa,Bq,B0 



uations can be characterized by few parameters and it is possible to solve those 
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For t, t' — > 00, jr << 1 and r = t — t', time translational invariance (TTI) holds 

lim C k (t,t) = q k , 

t—>-oo 

lim C k (t + T,t) = C k {r), 

t — >oo 

lim Cfe(r) = q k , 



(69) 
(70) 
(71) 



and 



lim R k (t + T,t) = R k (r). 



(72) 



In addition to the TTI regime, there is another long time non trivial regime, characterized by t, t' — * 00, fixing 
A = h{f) /h{t) and < A < 1, where the function hit) is an increasing function of t which the asymptotic analysis 
performed here is not able to determine. In this aging regime one has 







lim C k (t,h- l (Xh(t))) = q k C k (X), (73) 



t — *oo 



and 



lim q k C k {X) = qo,k, (74) 

A — »0 

lim C fe (A) = 1, (75) 



lim R k (t,Xt) = 7^ fc (A). (76) 
Also, for future convenience, it is useful to introduce the following order parameters: 

/dk 
7T 9k, (77) 

/dfc 
7T - 9fe ; (78) 

7T- Qo.k, (79) 

Z7T 

together with 

b = 2(q-q) , b Q = 2(q-q ). (80) 



X. EQUATIONS RELATING ASYMPTOTIC VALUES OF CORRELATION AND RESPONSE 

FUNCTIONS 



Using the ansatz discussed in section IX one can derive the following equations for C k (t,t') in the TTI regime: 

_,^dC k (T>) 



[T(p + fc 2 ) + d/dr] C fe (r) = 2TR k (-r) + ^V(b) \C k (r) - q k ] - | j* dr'V(B(r - r'))- 



dr' 

+2 f d P V{B{p))R k (p)+A [ d P V"(B(p))f(p)\c k (r) - q k C k (p) 
Jo Jo L 



(81) 



It is also possible to derive similar equations for R k (r) which, due to the Fluctuation Dissipation Theorem (FDT) 

1 dC k (r) 



R k {r) 



T d-, 



(82) 



are completely equivalent to Eq. (pl|). 

In the aging regime one gets the following equation for q k C(X): 



[T(p + k 2 ) -4 / dpV\B{p))f(p) 



q k C k (X) = 2 / d P V{B{p))R k (p) + -V'(B(X))(q k - q k ) 



T 



-4/ dpV"(B(p))f(p)q k C k (p/X)-4 dpV"(B(p))r(p)q k C k (X/p). 



(83) 



For R k (X) we obtain, 



T(p + k')-A I dpV"(B{p))f{p) 
o 



R k (X) = 

Again, one can see that both Eq. (|83|) and Eq. (fe4) can be solved by the ansatz 



l 

A P 



V"(B(X))r(X)(q k -q k )-4 I (B(p))f(p)R k (X/ p). (84) 
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(85) 



Eq. ( p5| ) is commonly referred to as a generalized FDT (GFDT). In principle, Eq. ( |S5| ) could have been written as 

4(a) = x -^m qk d ^ (86) 

1 d\ 

which could be applied to a many-step RSB scheme. However, as previously discussed, the present random het- 
cropolymer model can be identified with the random manifold problem with short range potential correlations. As 
such, it has one step RSB, and it is sufficient to use the simpler ansatz given in Eq. (B5|). 



For t — t' and t — > oo Eq.(^) gives 

T( f l + k 2 )q k = T+^V(b)(q k -q k ) + 2 j <I/>V'(B(/>) )/?/,(/•>) + 4 / ,/,,V"l 13 i ,, ) )rt ,,) [a, - y/l O. (//)]• IS7) 
Eq. ( |Sl| ) for t — > oo and then r — > oo results in 

T{fi + k 2 )q k = |v'(6)(& -q k ) + 2 f flfiVlBl/)) )/>',, (,>) i- 4 / tli>V"iB{p))r{p'hi), [I O I,,)]. | s.s) 



Also, Eq.(|8J) for A -> gives 

T(p + k 2 )q , k = 2V'(b ) f d P R k (p) + ^V(b )(q k -q k ). (89) 

Eqs. (|87|), (|88| ) and ( [89]) (and, equivalently, |8l], |8^ and |8~i|) contain TTI and aging parts. Thus, in principle, there 
are two anstze for solving them, leading to two phases: an ergodic one (without aging) and a glassy one (with aging). 



XI. ERGODIC PHASE 



Technically, assuming that aging is absent amounts to setting R k (X) — and C k (X) = 1 in (|87|), ( p8[ ) and (|8^ 

eading to the same 



(Equivalently, one could start from ( p7| ) and ( p8[ ) and exclude the aging part from the beginning, 
equations.) Thus, in the ergodic phase, equations (|S7|), (88) and (p9]) reduce to 



T(/i + k 2 )q k =T+ -V'(b)(q k - q k ) 

T(p + k 2 )q k = ^V'(b)(q k -q k ), 
2 

T(fi + fc 2 )(?o,fc = j, V'{b )(q k - q k ). 



(90) 
(91) 
(92) 



Note that ( |9l| ) and (|92|) enforce q k = q , k which is just equivalent to C k (X) = 1, so one gets only two equations. 
Solving them for q k and q k gives 



1 



qk - qk = 

qk 



p + k 2 
1 2 



V'(b)- 



After integrating over k and using 



dk 1 



p + k 2 T 2 y J (p + k 2 ) 2 ' 



dk 



2irp + k 2 2^/Jl ' J 2tt (fi + k 2 ) 2 4/z 3 / 2 ' 



(93) 
(94) 

(95) 
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we obtain 



2^/ 2 T 2 



q = 



2^J1 2/i 3 / 2 T 2 



v'(i/Vm). 



(96) 
(97) 



For T very small q and g blow up since the confinement term /xa;(s, i) 2 term becomes ineffective (see Eq. Q). For very 
large temperature q approaches zero but is never exactly equal to zero. 



XII. SPIN GLASS PHASE 



Keeping the aging parts and using the GFDT, Eqs. (f 



and (B9L) can be transformed into 



T{n + k 2 )q k =T+ -V(b)(l - x)(q k - q k ) + -V'(b )x(q k - q 0<k ) 



T 



T(n + k 2 )q k = -(V'(b)-xV(b ))(q k -q k ) + -V (b )x(q k - q 0<k ), 



T 



T 



T(» + k 2 )q , k = -V'(6 )(l - x)(q k - q k ) + -V'(b )x(q k - q 0>k ). 



Solving Eqs. @, @ and ( |10C| ) for q k , q k and q 0tk gives 

1 



qk - qk 



f i - 



k 2 + E ' 



1 1 



1 - x 



1 



qk - qo,k = - 



x /i + k 2 x // + &+£' 
2 1 
q k = (qk - qo,k) + ^2 V '(M + fc 2)2 > 



where 



X = x^(V'(b)-V(b )). 



Integration over k and using (051) gives 



b = 



1 



11 1 - x 1 

X y/Jl X ^//X + S' 

? = *> + 2^V(6o). 



Furthermore, Eq. (|3) with A = 1 gives 



i? fc (l)(M + fc 2 + S) = -(ft -ft) 



4V"(6) A 
r 

J>2 



and after using Eq. (101), integrating over k and using /j + £ = b 2 (see Eq. |l05| ) one gets 

= f(l) [T 2 + 6 3 V"(&)] . 



(98) 
(99) 
(100) 

(101) 
(102) 
(103) 

(104) 

(105) 
(106) 
(107) 

(108) 
(109) 



Eq. (|109|) with r(l) ^= implies marginal stability condition 
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T 2 = b 3 V"(b). 



(110) 



Also, equations (105) and (106) can be rewritten as 

V'(b)- V'(b ) _T 2 ^l(l \ 



&o 



6 

x 



2 fo V & 
6 ' 



(111) 
(112) 




lll| ) and ( |112| ) fully solve the model: ( |1 1C| ) gives 6 as function of T, (111) determines 6o as function of T 
determines x(T,fi) and Eq. (107) gives q(T,fi). Knowing b( T), b p(T, u) and q(T,u) determines q(T,/j,) 



and qo(T,u). Were we to impose the spherical constraint q = const, Eq. 107 could be used to relate fj, and T, and 
all order parameters could be expressed as functions of T only (q being fixed) . However, in this study we work with 
fixed T and \i allowing q to change. 



XIII. SOLVING THE EQUATIONS (PHASE DIAGRAM) 



The procedure of solving equations similar to the ones given in Eqs. ( |110| ), (111) and (|112| ) has been discussed in 
ref. p5fl . We apply a similar analysis to the random heteropolymer problem. In principle, three are three critical lines 
in the T, /i plane separating them (as shown in figure 1). 

Critical line (1): T = T max is the uppermost critical line (denoted in figure 1 by (1)); a bove this line Eq. flllGj ) 
has no solution. The value of T max can be determined from the graphical solution of Eq. ( |llOD depicted in figure 
2. Once T has been chosen (horizontal line labeled (T/T max ) 2 ) b is found from the interc ept o f (T/T max ) 2 line with 
— 6 3 V"(6) curve. From figure 2 it is clear that at b = b max the right hand side of Eq. (110) reaches a maximum; 



requiring ± [b 3 V"(b)] = gives 3V"(6) + bV"'(b) = and b max = Accordingly, T max = [-b 3 max V" (b max )] 



1/2 



Also, note that for fixed T, Eq. (110) has two solutions for b (denoted by 6/ and bu in figure 2). The first, physical 
solution (bi — > for T — > 0), in the interval [0, b max ] and second, unphysical solution (bu — * oo for T — ► 0), in the 
interval [b max ,oo). Accordingly, a model with a = (i.e. V(Ax) — 6 (Ax)) leads to an unphysical phase diagram, 
since for a — » physical branch [0, b max ] disappears (b max —* 0). 

Clearly, the form of V(b) for small b has to be modeled carefully and the choice V(Ax) = S(Ax) simply fails in 
that respect, giving V(0) = oo. Thus when formulating the problem, if there is to be a possibility of freezing at low 
temperatures (b — > as T — > 0), the bead-bead interaction V(Ax) has to be regular for small Ax. Similar small 
distance regularization problem of bead-bead interaction was encountered in ref |2lj . 

Critical line (2): corresponds to b = bo. Fro m Eq . (112) it follows that b — bo = l/y/JI. The equation of the 
critical line is obtained by inserting b = 1/^/71 into ( |110| ): 



(T/Bf 



V 3/2 



-1/2 



-(7 + 1) 



(113) 



U E [l-i r , 



where \i max solves Eq. (113) with T = T, 



The value of x c at the critical line can not be directly obtained from Eq. (112). Instead, one has to approach the 



critical line and obtain the lim iting value of x: for example, first, one assumes that point (T c ,[i c ) is at the critical 
line (T c and fi c satisfy Eq. |l!3| ) and, then, T(e), /i(e), 6(e), bo(e), x(e) approa ch th ei r val ues at the critical line for 
e — * 0. Naturally, the dependence on e has to be chosen consistently with Eqs. (|110|), (|1 1 1|) and Qll2|). Since one has 



five variables and three equations which relate them, two variables have to be specified as, e.g., bo(e) = b c + e, with 
be = 1/ y/ Uc., and T(e) = T c . The other three variables 6(e), /i(e) and x(e) have to be determined from (11C), (111) 
and ([ill ): 



V'(b c )-V(b c + e) _ T 2 ^/JI(e) (I 



(114) 
(115) 
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Equation (110) is trivially satisfied and does not enter the analysis. At first order in e Eqs. (|115|) and (114) give 



x(0) 



2/i! 



1 ,, n , ,, n v 2V W (6 C ) 
rA* (0) , M (0) - 



3 I?v^ 



3/2' 



(116) 



which, together with T c 2 = — /i c 3 ^ 2 V"(6 C ), gives the value for x at the critical line (2), 

_ 1 V"'(b c ) 
Xc ~ 3V"(6 C )' 



(117) 



Using the explicit form for V gives 



7 + 1 b c 
3 b c + a 



(118) 



with b c on the critical line. b c — b max gives x c = 1 while for fe c = one gets x c = 0. 

Thus, at the critical line (2), close to T max , x c is very close to 1 and as T (fi) decreases (increases) x c drops to zero. 
Also, at the critical line (2), the transition to the ergodic phase is continuous in b and bo and discontinuous in x. 

C riti cal line (3): at this line x = 1 and Eq. flll2j ) gives b = 1/JJi > b. The equation for this critical line is given 

by ra. 



V(b) - vwyg _ J 12 

1/y/Jl-b 2 6 



(119) 



Once T is chosen, 6 is determined from (110) and upon solving Eq 



(3) is depicted in figure 1, where it was obtained by solving Eq. (119) numerically. The line starts from (/x 

la 



and then drops to (0,T*) where T* is given from Eq. (110) with b = b* and b 



L 1 9 ) one obtains fi as function of T. Critical line 

) 
as 



can easily be checked by inserting those assumptions in Eq. (119). Also, bo 



7-2- 

b as /i — 



T 

max : J max ) 

Thus, as b — > b* , bo — > oo, 
Mmaz- Thus, contrary to (2), 



on line (3) the transition to the ergodic phase is discontinuous in b and bo while continuous in x. 

Also, for arbitrary /i, when T gets close to zero b approaches and bo grows to infinity. This simply means that for 
low temperatures the heteropolymer freezes completely: x(s,t + t) = x(s,t) for arbitrary r and t sufficiently large. 
On the other hand, for fixed T and vanishing /i, Eq. (Ill) gives bo — + oo, while b stays fixed by Eq. (11C). 

For small fi Eqs. ( j96"U9l| ) give q/q oc ^( d ~ 2 )/ 4 . Thus, for fi — one gets q/q — 0. Also, as discussed in the preceding 
paragraph, in the glass phase for ^ — ► one has bo — + oo and b — const, which gives q/q = 1. Thus contrary to the 
ergodic phase, where vanishing fi lead to paramagnetic-like behavior, in the glass phase the system gets trapped in 
one of many states separated by diverging barriers. Interestingly enough, adjusting /j, so that the radius of gyration 
R g scales according to Rg ~ N and using the relation R 2 g ~ l/^//I (which is exact for the Gaussian coil) gives 
/i oc N~ 4 / d and q/q oc N~( d ~ 2 ^ d . Thus, in the thermodynamic limit q/q — > 0. 



XIV. DISCUSSION 



We have presented a detailed derivation of the equations of motion of a random heteropolymer using SUSY formalism 
and a Gaussian variational ansatz. In deriving these equations, we have used a long-chain approximation, considerably 
simplifying the dynamical action. Furthermore, imposing translational invariance, we have shown that, as happens 
in statics, within the Gaussian variational ansatz the equations of motion become identical to those for a manifold in 
a random potential with power law correlations. 

Clearly, this result is strongly related to the particular variational ansatz employed here, and its generality beyond 
this framework remains an open question. Nevertheless, the existence of this mapping at the level of the GVA is rather 
intriguing. It connects the random heteropolymer model with many physical systems, such as a manifold pinned by 
impurities, interfaces in a random field, the glassy phase of vortices in high-T c superconductors, directed polymers 
in a random potential, and surface growth on disordered substrates. It would be interesting to understand to what 
extent the mappings to these problems extends beyond the GVA. 
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By making the standards 1RSB aging ansatz for response and correlation functions we found the asymptotic solution 
of the dynamical equation. The validity of this ansatz has been carefully checked elsewhere: in the context of random 
manifold problem it was shown that one step replica symmetry breaking ansatz can be used to describe random 
manifold with short range correlations, and we have applied this results to the random heteropolymer. 

The analytic solutions show that, as expected, the random heteropolymer has characteristic properties of spin glass 
systems: aging and ergodicity breaking. Furthermore, the dynamical phase diagram is different from that for statics. 
In dynamics starting from a random condition, the polymer get stuck at energies higher then the ones of the native 
state. 

In a more realistic approach to heteropolymers, we expect that finite dimensional, and finite length chain effects 
will be responsible for ultimate restoring of ergodicity. Our study should be taken as an indication of a time regime 
where trapping effect and aging could be observed. 

One of the motivations for this paper, mentioned at the beginning of the introduction, was the hope that it might 
provide some insight into the dynamics of proteins, including their folding. However, it is fairly well understood by 
now that protein dynamics are influenced strongly by the existence of an energetically favored native state, a feature 
absent from the random heteropolymer model we have studied here. In work currently in progress, we are extending 
the analysis presented here to models in which the two-body interactions -B s . s ' are systematically biased, with a 
tunable strength, to favor particular "native" states. Such models provide an opportunity to study the competition 
between the attraction to a native state and the glassiness produced by the randomness and frustration. (Refs. (3?] ficj ] 
treat equilibrium aspects of this competition.) 
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FIG. 1. Phase diagram of dynamic random heteropolymer model in fi, T plane. Critical lines are denoted by (1) T = T ma x; 
(2) b = bo = /i -1 / 2 {x ranges from 1 to 0); (3) x = 1, bo = /jT 1 ^ 2 > b. Below (2) and (3) lies glassy phase and above ergodic 
phase 
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FIG. 2. Graphical solution of Eq. ( |l 10[ ) ■ Equation has two solutions for T < Tmax denoted by "I" and "II". For T = T max 
there is only one solution b — bbmax- Solution I is physical and solution II is unphysical. 
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